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A thermal- model analysis of particle production of p — p collisions at ^/s — 17 GeV using the 
latest available data is presented. The sensitivity of model parameters on data selections and model 
assumptions is studied. The system-size dependence of thermal parameters and recent differences in 
the statistical model analysis of p — p collisions at the super proton synchrotron (SPS) are discussed. 
It is shown that the temperature and strangeness undersaturation factor depend strongly on kaon 
yields which at present are still not well known experimentally. It is conclude, that within the 
presently available data at the SPS it is rather unlikely that the temperature in p — p collisions 
exceeds significantly that expected in central collisions of heavy ions at the same energy. 
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I. INTRODUCTION 

The statistical model has been used to describe par- 
ticle production in high-energy colhsions for more than 
half a century [l[ . In this period it has evolved into a very 
useful and successful model describing a large variety of 
data, in particular, hadron yields in central heavy-ion 
collisions 0, 0] have been described in a very systematic 
and appealing way unmatched by any other model. It 
has also provided a very useful framework for the cen- 
trality 0] and system-size dependence 0, Q of particle 
production. The applicability of the model in small sys- 
tems like p — p 01 and e"*" — e~ annihilation @ has been 
the subject of several recent publications [ol. [lOl. [lH . 

The statistical-model analysis of elementary particle 
interactions can be summarized by the statement that 
the thermal parameters show almost no energy depen- 
dence in the range of -/s =14 - 900 GeV with the tem- 
perature being about 165 MeV and the strangeness un- 
dersaturation factor 75 being in the range between 0.5 
and 0.7. 

In the context of the system-size dependence of particle 
production, the p—p collisions at ^/s — 17 GeV have been 
analyzed in detail recently. Based on similar data sets, 
the extracted parameters in different publications devi- 
ated significantly from each other: in a previous anal- 
ysis (Rcf. [1, El) we derived T = 164 ± 9 MeV and 
7S = 0.67 ± 0.07, with x^/n 1-7/3, while the authors 
in Ref. @ obtained T = 178 ± 6 MeV, js = 0.45 ± 0.02 
with x^/" = 11/7. These findings motivated different 
conclusions: In Ref. Q no system size dependence of the 
thermal parameters was found, except for 7s which tends 
to increase when more nucleons participe in the collisions 
but this rise is weaker than the errors on the strangeness 
suppression parameter. In Ref. it was therefore con- 
cluded that the hadron gas produced in central collisions 



at -^/i = 17 GeV reaches its limiting temperature. Based 



on Ref. [6| on the other hand, it was argued in Ref. [12 1 
that, decreasing 75 and, in particular, increasing tem- 
perature towards smaller systems allow for probing QCD 
matter beyond the freeze-out curve established in Pb-Pb 
and Au-Au coUisions [l3l.[T4j. 

The goal of this paper is to understand the origin of 
these rather different thermal-model results obtained in 
the analysis ol p — p data. We use an up-to-date com- 
plete set of data and discuss the sensitivity of the thermal 
model parameters on their values. We present systematic 
studies of data used as inputs and the methods applied 
in their thermal model analysis. 

The paper is organised as follows: In Section |TT] we 
discuss the experimental data on which different analysis 
are based. In Section Hill we summarize the main features 
of the statistical model and present the analysis of the 
SPS data obtained in p — p collisions. In the final section 
we present our conclusions and summarize our results. 



II. DATA 

The data used throughout this paper for hadron yields 
in p — p collisions at -^5=17. 3 GeV are summarized in 
Table HI Data in column Set A were exploit in our pre- 
vious analysis [Tlj and the corresponding references are 
given in the table. If the numerical values deviate in the 
analysis of Ref. @ , they are listed in column Set B. The 
relative differences between the particle yields from sets 
A and B are also indicated in Table [D The commonly 
used data in the statistical model description of particles 
production in p — p collisions at the SPS are displayed 
below the horizontal line in Table HI 

In Table In] the experimental data are grouped in sets 
which are used in Section IIIII to perform the Statistical 
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FIG. 1; Charged kaon yields (left panels), negatively charged 
hadron h~ and pions n~ (upper right) and K~/tt~ ratios 
(lower right panel) in p — p collisions as a function of labo- 
ratory momentum. The charged kaons, K'^ (diamonds) and 
K~ (crosses) yields are from Ref. [S^]. The lines are fits to 
data. The SPS yields from Ref ^ (circle) and from Ref. ^ 
(triangle) are also shown. The negatively charged hadrons are 
from Ref jH). 



TABLE I: Particle yields (47r integrated) in minimum bias 
p — p collisions at ^/s = 17.3 GeV. Numerical values of Set A 
are from Ref. [TJ] . For Set B and common values (data below 
the horizontal line) the references are given in the last column. 



Particle 


Set A 


Set B 






Ref 


7V + 


3.02 ± 0.15 


3.15 ± 0.16 


4.4% 


10.5% [15] 


■K~ 


2.36 ± 0.11 


2.45 ± 0.12 


3.8% 


5.9% 


[15] 


K+ 


0.258 ± 0.055 








[16] 


K- 


0.160 ± 0.050 








[16] 


K+ 




0.210 ± 0.021 


19% 


62% 


[18] 


K- 




0.130 ± 0.013 


19% 


74% 


[18] 


A 


0.116 ± 0.011 


0.115 ± 0.012 


0.9% 


9.1% 


[la 


A 


0.0137 ± 0.0007 


0.0148 ± 0.0019 


8.0% 


171% 


[12 


Kl 


0.18 ± 0.04 






[18] 


P 


0.0400 ± 0.0068 






[18] 


A* 


0.012 ± 0.003 






[19] 


i 


0.0120 ± 0.0015 






[20] 




0.0031 ± 0.0003 










0.00092 ± 0.00009 










0.00026 ± 0.00013 






m 


n+ 


0.00016 ± 0.00009 






m 



Model analysis. In the following we motivate the partic- 
ular choice of data in these sets and discuss how they can 
influence the model predictions on thermal conditions in 
p — p collisions. 

The data set Al is most restricted. Firstly, the produc- 
tion yields of S and 17 are not included because their nu- 
merical values are only preliminary (Ref. [2]|). Secondly, 



TABLE II: Different sets of particle yields used in the thermal 
model fits. The type A sets contain numerical values from 
the left (and common) columns of Table |T] The type B sets 
contain data from Set B and common columns of Table |l] 



Set Particles 








Comment 


Al 


^± 


K± 


A A p 








Set of Ref [11] 


A2 


^± 


K± 


A A Kg p 


A* H- H+ 




n+ 


A3 


^± 


K± 


A A Kg p 


A* <t> 








A4 


^± 


K± 


A A Kg p 








K± from Ref [181 


Bl 


^± 


K± 


A A Kg p 


A* 






A* contribution 


B2 


^± 


K± 


A A Kg p 


A* 


Q.- 


n+ 


Fit A of Ref. [6] 


B3 


^± 


K± 


A A Kg p 


A* 






Fit B of Ref. [6| 


B4 


^± 


K± 


A A Kg p 








Set Bl without A* 



the A* resonance is also not included so as to restrict 
the analysis to stable hadrons. Finally, the cj) meson is 
omitted in Set Al since this particle is difhcult to address 
in the statistical model due to its hidden strangeness as 
discussed in Ref. @. 

The lower yields of charged kaons in Set B of TableUare 
taken from results published in conference proceedings 
[l^. Such kaon yields are in disagreement with trends 
from data measured at lower and higher energies as seen 
in Fig. m 

The left panels of Fig. [T] show the charged kaon multi- 
plicities from p — p interactions at lower and higher beam 
momenta together with data from Table HI The lines 
in this figure are simple parametrizations interpolating to 
SPS energies. The K~ yield from [l^ is seen to be 7% 
below the expected value from the above paramctriza- 
tion, however agrees within errors. The abundance 
from [l^ is by 24% lower and its error is only 10%. As 
we discuss below, such a low value for the multiplicity 
of charged kaons influences the statistical model fit in an 
essential way. 

The upper right panel of Fig. [1] shows the negatively 
charged hadrons froni p — p interactions at several beam 
momenta from Ref. [2^. As indicated in Ref. [1^ the 
K~ , p and E~ supplement the 7r~ yield. In this case, the 
ratio Tr~ /h~ amounts to 91%. Including more sources 
of feed-down, the tt~ /h~ ratio stays at the same level 
as long as A and Kg can be separated. Consequently, 
to calculate the negatively charged pions from yields 
one can use the above 91% scaling factor. Figure [T] (top 
right panel) shows the fit to yields as a function of 
beam momenta and then by rescaling the expected result 
for the puh dependence of the negatively charged pions. 
The yields of tt~ at SPS from Table [J agree quite well 
with that expected from an interpolation line shown in 
Fig. [1] They are only slightly higher, by 1% for yields 
taken from [l^] and by 5% for yields used in Ref. [y]. 

The lower right panel in Fig. [1] shows the K^/tt" ra- 
tio at SPS compared to the interpolated data from other 
beam momenta. The mean value of the K~ /Tr~ used in 
[m is 8% below the interpolated line but agrees within 
errors, while the corresponding value used in is 28% 
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smaller and exhibits an error of only 11%. Clearly, the 
above differences in the K~ /tt^ ratios influence the ther- 
mal model fits. 

In general, a smaller kaon yield implies a stronger sup- 
pression of the strange-particle phase space resulting in 
a smaller value for the strangeness undersaturation fac- 
tor 7s. If other strange particles are included, then the 
strong suppression caused by 7s has to be compensated 
by a higher temperature. This might be one of the origins 
for the different thermal fit parameters obtained in Rcfs. 
[m and . In order to quantify this we have selected a 
data set A4 which is equivalent to the Set Al but with 
the kaon yields of Ref. being replaced by the values 
from Ref. [H. 

The Set Bl is (besides the A*) equivalent to Al but 
with numerical values for particle yields from column B in 
Table HI The Set B4 is used to demonstrate the influence 
of the A* resonance on thermal fit parameters. The Sets 
A3, B3 and A2, B2 are chosen to study the influence of 
the (j) meson and the multistrange hyperons on thermal 
fit parameters. 



III. STATISTICAL MODEL ANALYSIS 

The usual form of the statistical model formulated in 
the grand-canonical ensemble cannot be used when either 
the temperature or the volume or both are small. As as a 
rule of thumb one needs VT^ > 1 for a grand-canonical 
description to hold [23, HI]. Furthermore, even if this 
condition is matched but the abundance of a subset of 
particles carrying a conserved charge is small, the canon- 
ical suppression still appears even though the grand- 
canonical description is valid for the bulk of the produced 
hadrons. There exists a vast literature on the subject of 
canonical s upp ression and we refer to several articles (see 
e.g. iSHl). 

The effect of canonical suppression inp—p collisions at 
ultra-relativistic energies is relevant for hadrons carrying 
strangeness. The larger the strangeness content of the 
particle, the stronger is the suppression of the hadron 
yield. This has been discussed in great detail in ^28j . 

In line with the previous statistical model studies of 
heavy-ion scattering at lower energies, the collisions of 
small ions at SPS revealed Q that the experimental data 
show stronger suppression of strange-particle yields than 
what was expected in the canonical model [1, [H, Isd^ . 
Consequently, an additional suppression effect had to be 
included in order to quantify the observed yields. Here we 
introduce the off-equilibrium factor 75 < 1 which reduces 
densities Ug of hadrons carrying strangeness s by — > 

ns-7s' lii- 

We investigate whether or not all quantum numbers 
have to be conserved exactly in p — p collisions within a 
canonical approach by comparing data with two model 
settings: 

• Canonical (C) Model: all conserved charges, i.e. 



TABLE III: A list of parameters needed to quantify particle 
yields in the strangeness canonical (SC) and canonical (C) 
statistical model (see text). The symbols S, Q and B are 
the strangeness, electric charge and baryon number respec- 
tively with for i = {S,Q,B) being chemical potentials re- 
lated with conservation of these quantum numbers. 7s is the 
strangeness undersaturation factor, T the temperature and 
R the radius describing the spherical volume of the collision 
zone. 





SC 


model 


C model 


Fit parameter 


MB 


R 


S Q R 


Constrained param. 


Mq: 


B/2Q = 0.5 




Fixed param. 






B = 2 


Fit / Scan param. 


T 7s 


T-ys 


No. of parameter 


5 




6 




Fit 


X'' scan 


Fit scan 


No. of free param. 


4 


2 


5 3 


No. of fixed param. 


1 


3 


1 3 



strangeness, electric charge and baryon number are 
conserved exactly within a canonical ensemble. 

• Strangeness Canonical (SC) Model: only 
strangeness is conserved exactly whereas the 
baryon number and electric charge are conserved 
on the average and their densities are controlled 
by the corresponding chemical potentials. 

The parameters of these models are listed in Table IIIII 
In the following we compare predictions of the above sta- 
tistical models with p — p data summarized in different 
sets discussed above. 



A. Comparative study of p — p data at SPS 

We start from the analysis of data set Al and mod- 
ify it stepwise to find out in which way one matches the 
conclusion of larger temperature in p — p than in central 
A — A collisions at SPS as indicated in Ref. All nu- 
merical values of model parameters are listed in Table [TVl 
A detailed discussion on their choice and correlations is 
presented in the Appendix based on the x^/n systcmat- 
ics. 

The fit to data set Al in the SC model complies with 
our previous analysis from Ref. [ll[, see also Fig. [2] (top). 
The SC model fit to these data does not change when 
including A* hyperons resulting in the same values of 
thermal parameters and their errors as summarized in 
Table HVl 

The most striking effect on thermal parameters is ex- 
pected when replacing the kaon yields in Set Al (Fig. [31 
top) by those from Ref. [Hi, Set A4, Fig. [U (bottom). 
Indeed, smaller kaon yields cause an increase of temper- 
ature and a decrease of 7s. These changes come along 
with a reduced volume and in case of the SC fit with in- 
crease of the baryon chemical potential. The kaons from 
Ref. [13 dominate the fit because their errors are 10% 
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TABLE IV: Thermal parameters extracted within the 
strangeness canonical (SC) and canonical (C) model (see 
text) from fits to 47r-integrated data in p — p collisions at 
-^5 — 17.3 GeV. In the SC model analysis for data set Bl 
a fit does not converge, the minimum of the scan is dis- 
played. The fits to data sets A3 and B3 do not exhibit a 
minimum in the parameter range considered, thus only the 
tentative parameters of a possible minimum are indicated. 



SC model results 



Set 


r(MeV) 


75 


R (fm) 


fiB{MeV) 




Al 


163 ± 5 


0.68 ± 0.05 


1.50 ± 0.11 


208 ± 14 


1.7/4 


A2 


168 ± 1 


0.66 ± 0.02 


1.37 ± 0.03 


221 ± 8 


8.6/9 


A3 


> 190 


^ 0.5 


< 1.1 


> 250 




A4 


177 ± 5 


0.59 ± 0.03 


1.23 ± 0.10 


233 ± 16 


5.1/4 


Bl 


176 


0.56 


1.24 ± 0.01 


240 ± 12 


7.7/7 


B2 


179 ± 5 


0.61 ± 0.02 


1.19 ± 0.09 


242 ± 18 


16/9 


B3 


> 190 


^ 0.5 


< 1.1 


> 250 




C model results 


Set 


r(MeV) 


7s 


R (fm) 




x/n 


Al 


175 ± 5 


0.57 ± 0.04 


1.33 ± 0.09 




0.5/3 


A2 


174 ± 4 


0.59 ± 0.02 


1.34 ± 0.08 




6.6/8 


A3 


189 ± 5 


0.46 ± 0.02 


1.12 ± 0.09 




23/5 


A4 


181 ± 4 


0.52 ± 0.03 


1.22 ± 0.07 




3.5/3 


Bl 


177 ± 5 


0.51 ± 0.03 


1.30 ± 0.09 




6.8/4 


B2 


180 ± 4 


0.56 ± 0.02 


1.23 ± 0.08 




18/8 


B3 


178 ± 5 


0.45 ± 0.02 


1.30 ± 0.10 




19/5 


B4 


177 ± 5 


0.50 ± 0.03 


1.31 ± 0.09 




4.8/3 



while the uncertainties of the and K~ yields, taken 
from Ref. [l^, are 21% and 31%, respectively. Conse- 
quently, the smaller errors dominate the statistical model 
fit. 

In the next step we add A*, S and Q hyperons result- 
ing in Set A2, see Fig. |31 middle panel. The measured 
hyperon multiplicities coincide with the model results ob- 
tained before, thus within errors the statistical model 
parameters remain unchanged. We focus on Set A3 and 
add the cf) meson. Fig. [3] bottom. In this case the tem- 
perature is indeed much higher. In the SC model the 
thermal parameters obtained from the Sets A3 and B3 
appear to be meaningless and unphysical due to the cf) 
meson contribution. 

Additional to different kaon data in Sets A and B there 
are also slightly different values for pions and A yields, 
see Table ID Compared to results obtained from data 
Sets A4, in the fit of Set Bl the higher pion yield re- 
duces the strange to non-strange particle ratios resulting 
in slightly smaller value of 75. The fit obtained with the 
Set B4 yields similar results as that obtained from Set 
Bl indicating that the A* resonance docs not affect the 
model parameters. 

In the canonical (C) model analysis the Sets Al and A2 
tend towards a slightly higher temperature and smaller 
7S than that obtained in the SC analysis. The situation 
is different for Sets A3 and B3 that include the (j> meson. 
Here in the SC model the temperature is very high and 
7S ~ 0.5. In the C model the temperature decreases and 




0.4 0.5 0.6 0.7 

's 

FIG. 2: The scan in the (T-7s)-plane. Starting from its 
minimum, increases by 2 for each contour line. Upper 
figure: fit to data set Al in the model where only strangeness 
is conserved exactly (SC). Lower figure: fit to data set A4 in 
the canonical (C) model. The minima are indicated by the 
crosses. 



7s drops below 0.5. We can conclude that in the case 
where the (f> meson is included in the fit, one needs to 
apply the C analysis to get lower temperatures, however 
with very small values of 75 and a large x^ jn. For Set B3 
the numerical results for T and 7^ summarized in Table 
IV coincide with that obtained in Ref jH, however with 
a larger x^ jn [Slj . 



IV. DISCUSSION AND SUMMARY 

The statistical-model analysis of hadron yields for 
•p — V colHsions at ^/s = 17 GeV from Refs. Q and 0], 
yield different results and lead to different conclusions on 
the system-size dependence of thermal parameters [5l.[l^. 
In this paper we have reanalyzed the p — p data and stud- 
ied the sensitivity of the thermal fit to data selection and 
on model assumptions. We have shown that different 
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FIG. 3: As in Fig. [2] but for data Set Al (top), Set A2 
(middle), Set A3 (bottom). Canonical ensemble. 

conclusions from Refs. Q and Q are mostly due to dif- 
ferences in data selections. 

Slightly different numerical values for charged pions 
and A hyperons used in Refs. Q and @ as well as the 
contribution of the A* resonance altered thermal param- 
eters only within errors. However, the used charged kaon 
yields in both approaches differ substantially. We have 
argued that data of kaon yields in Ref. 0] deviate from 



trends seen in data at different energies resulting in a 
higher temperature. 

We have shown that higher kaon yields expected from 
the systematics in the energy dependence in p — p col- 
lisions are in line with data on multi-strange baryons. 
Unlike the hyperons. when adding the (f> meson the ther- 
mal model fit leaves a reasonable range of parameters 
resulting in a very high temperature exceeding 190 MeV 
and large ^ jn. We have quantified the modifications 
of these results when including an exact conservation of 
all quantum numbers in the canonical statistical model. 
We have shown that in the absence of meson the ther- 
mal fits are rather weakly influenced by canonical effects 
due to an exact conservation of the baryon number and 
an electric charge leading in some cases to a system- 
atic increase of the freezeout temperature. Fits including 
the (/) meson are sensitive to an exact conservation of all 
quantum numbers resulting in lower temperatures. How- 
ever, the thermal model analysis of data sets with hidden 
strangeness has the largest X^/n indicating that this par- 
ticle cannot be addressed properly in this model. 

From our analysis, we conclude that within the 
presently available data an p — p collisions at SPS en- 
ergy and uncertainties on thermal parameters obtained 
from fits within the statistical model, it is rather unlikely 
that the temperature va. p — p collisions exceeds signifi- 
cantly that expected in central collisions of heavy ions at 
the same energy. 
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APPENDIX A: THE CONTOURS OF 
THERMAL FITS 

In this appendix we quantify the choice of thermal 
parameters within the statistical model through x^- 
contours in the parameter space. Since the temperature 
T and 75 arc of particular interest here the quality of the 
fits are shown in Figs. 2 and 3 in the (T — 75)-planc. In 
these figures for fixed (T, 75)-pair the remaining model 
parameters were fitted and the corresponding ^ was cal- 
culated. 

Figure [5] (top) shows the analysis of the data Set Al 
within the strangeness canonical (SC) model. The analy- 
sis in the model with canonical treatment of all conserved 
charges (C) is shown in Fig. [3] for all data sets besides 
the Set A4 which is presented in Fig. [2] (bottom). 

In the C model description of data Set Al there is a 
large region of a very low ^ which manifests the expected 
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anti-correlation of T and 75. Reasonable fits are possible 
over a large range of parameters. For the Set A2 the 
minimum is located at the same temperature and slightly 
higher 75 . The contributions of S and Q baryons disfavor 
small values of 75. 



The (j), as seen in Fig. 3, directs fits towards very 
high temperatures and very strong strangeness suppres- 
sion. Also the pattern of (T-^s) anti-correlations shows 
decreasing with increasing temperature at fixed 75. 
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